Coulomb displacement energies, energy differences and neutron skins. 
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A Fock space representation of the monopole part of the Coulomb potential is presented. Quan- 
tum effects show through a small orbital term in l{l + 1). Once it is averaged out, the classical 
electrostatic energy emerges as an essentially exact expression, which makes it possible to eliminate 
the Nolen-Schiffer anomaly, and to estimate neutron skins and the evolution of radii along yrast 
states of mirror nuclei. The energy differences of the latter are quantitatively reproduced by the 
monopole term and a schematic multipole one. 
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The electrostatic energy of a sphere of radius R and 
charge Ze is easily calculated to be Ecs = 3e^Z^/5R 
{Ecs stands for "simple Coulomb energy" ) . It is un- 
der this guise that the Coulomb field enters the Bethe- 
Weizacker mass formula, and becomes a basic quantity 
in nuclear structure. Direct evidence of entirely Coulomb 
effects has long been available from displacement ener- 
gies betvifeen mirror (or in general, analog) ground states 
(CDE) and more recently from differences in yrast 

excitation energies in mirror p/-shell nuclei (CED) 

The CDE range from few to tens of MeV. They 
should be given mainly by Ecs, but are underesti- 
mated by mean-field calculations: the Nolen-Schiffer 
(NS) anomaly still an open problem. 

The CED are very small (of the order of 10-100 keV). 
They have been semiquantitatively explained by shell 
model calculations using Coulomb matrix elements but 
there is some room for improvement. 

The present letter intends to give a unified microscopic 
description of CDE and CED, by separating the Coulomb 
field into monopole and multipole components Vc = 
e^{Mu/hf/^lr = Fcm + Vcm following ref. i(DZ). 
The monopole Vcm contains all terms quadratic in scalar 
products of Fermion operators a\ ■ aj . Its diagonal part 
involves only proton number operators, and should be re- 
sponsible for Ecs, and hence for the observed CDE. The 
non diagonal part will not be considered here: it leads 
to isospin mixing, but energetically its effect is small. 
The multipole Vcm contains all non- monopole matrix 
elements and accounts for much of the CED. 

As an introduction we examine the origin of the NS 
anomaly, by calculating Ecs with proton radii of the form 

R. = ro (^1 _ _ + V. (1) 

{A~N + Z,t = N~Z,Tr,iy = protons, neutrons.) 

Isospin conservation is assumed, which implies that R^, 
is the same as i?^ for the mirror nucleus, obtained by in- 



terchanging N and Z. Therefore, in Eq. (|i|), for N > Z, 
V > represents a uniform contraction of the two fluids, 
while C > represents a 7r-contraction and a i/-dilation. 
Hence, a (fractional) neutron skin can be defined as 

^(0 ^ , J- f;,, ^ N>Z^ (2) 

The exponential factor takes care of the increase in Rt, 
observed in the light nuclei. 

2? is a phenomenological term that accounts for shell 
and deformation effects. It is a sum of two quar- 
tic forms XS^^Si, + plQttQi/ that vanish at the spin- 
orbit (or EI: extruder-intruder) closures at or Z 
6, 14, 28, 50, 82, 126. Defining = (p^-hl)(p^ + 2) + 2, 
the degeneracy of the EI^ shell (e. g., p^r = 3 for Z = 28 
to 50), DrT, = pTriPTT + 1), the degeneracy of the non- 
intruder subshells; the factors 8-r6 /S'tt — ziDT,-z)/Dl and 
— ^iEnr — z)/D'^ (z= number of valence protons). 
The parametrization is a variant of the ones in P,p^ . 

Fits to i?7r = ■\/5(r2) /3 — where (r^) is the measured 
mean square radius — for 634 nuclei (with N > Z except 
two cases!) yield (rg and rmsd in fm) 
ro=1.236 w=0.00 C=0.94 A=6.2 ^=14.6 rmsd=0.018: 
A good fit with a Huge Skin (HS hereafter). 
ro=1.220 u=0.61 C==0.00 A=5.7 /^=27.0 rmsd=0.012: 
A much better fit with Zero Skin (ZS hereafter). 
ro=1.226 u=0.45 C=0.29 A=5.7 Ai=24.0 rmsd^O.Oll: 
An even better fit with a Minute Skin (MS hereafter). 

Adding an exchange term to Ecs, replacing Z^ by 
Z{Z — 1) (conceptually better, as Vc is two-body), the 
main contribution to the monopole energy takes the form 

dZ{Z -l){l-cZ~^) 3e2 0.864 

Ecm — , « = -z — = , (3) 

P 5ro To 

d in MeV and p = (i?7r — T^)!^^: The V correction can be 
left out for simplicity as it does not affect what follows; 



1 



the high quahty fits become average ones (rmsd trebled), 
but still sufficient for our purpose. For the exchange term 
the usual choice is /i = 2/3, while c varies appreciably 
Here we set ft, = 1 (explained in and after Eq. (||)). The 
overall factor d is fixed. Nevertheless, it will be allowed to 
vary, under the name dj, together with c, in a fit to the 
183 available CDE from reference The values 

are calculated by adding 200 keV to the experimental 
errors, to account for uncertainties in the calculations. 
Obviously, consistency between and Ecm demands 
df = d. The three fits to Rt^ lead to: 

d=0.699 rf/=0.77 c = -0.5 w=0.00 C=0.94 x^=0.91: 
HS leads to a large overestimate of d. In other words: 
keeping the correct d leads to a large underestimate of 
the CDE: the NS anomaly 12|, unresolved to this day 
(see however |p^). 

d=0.708 d/=0.69 c= 1.3 u=0.61 C=0.00 x^=l-46: 
ZS leads to a small underestimate of d. 

d=0.705 d/=0.71 c= 0.9 w=0.45 C=0-29 x2=i.20: 
MS leads to d/ « d. The NS anomaly disappears. 
Clearly, ( ^ 0.3 is a good guess in Eq. (||). 

The NS anomaly occurs because mean field calcula- 
tions that yield good i?^ systematically overestimate R^: 
Ref. |lj] contains a nice illustration of the problem. It is 
ironical to note that a small neutron skin was recognized 
as a possible solution of the anomaly — but rejected — by 
Nolen and Schiffer |^|. Many reasons explain why this 
rejection held for so long. We retain only two: neglect 
of the v{t/A)'^ term in Eq. (0), and lack of confidence in 
the validity of Ecm- As we show next, Ecm in Eq. (|^) 
must be trusted, as it is basically an exact form of Vcm- 

By definition, the diagonal monopole part of Vc is 
([J] = 2J+1): 



i<k 



5ik) 
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Vk 



(4) 



The label k = plj stands for the quantum numbers spec- 
ifying a given harmonic oscillator (ho) orbit {p is the 
principal quantum number). Restricting the sum to the 
first K major shells containing r orbits, Vcm is brought 
to a sum of factorable terms by diagonalizing the matrix 
i{Vifc} through the unitary transformation U: 



Vcm = 



E 



(5) 



where C — '^^^ ZnV„n/2 is the one-body counterterm 
in Eq. (^) left out of the diagonalization. By rescaling 
£n — TEn, Ukn — {''')~^^^Ukn, the rcsults become inde- 
pendent of r. To fix ideas choose k = 8, i. e., r = 36. We 
expect to extract something close to the operator, in 
which case only the highest eigenvalue (E^q) should be 
non vanishing, with Ukse = 1- The diagonalization pro- 
duces indeed an E^q = 0.192 that is 30 times larger than 



the next and over 100 times larger than the second next. 
(Increasing the number of shells would only increase the 
number of negligible terms: Eq (^ is an exact represen- 
tation in Fock space.) Fig. |^ shows that {Uk = Uk^e ) is 
very well approximated by the form Up — 0.01 L^, 
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FIG. 1. Dominant factor in the DZ decomposition of Vcn 



where Up ~ J^ji'^j + '^)Uk is the average of Uk over j- 
values, and — l{l + 1) — p{p + 3)/2 is an l{l + 1) term 
referred to its centroid. This result is still almost exact. 
The presence of is interesting, but here we average it 
out, to transform the operator Zp Up into a c-number 
by taking expectation values in Eq. (|^) over ho closed 
shells. Introducing b = [h/MuY^'^ , we find 



{Vera) « y 



E:i,{Y,ZpUpf -{C) (6) 



^36 $^(P+1)(P + 2)C/J -(C) 



« ^0.445 {Z{Z - l))(i-i/i2)/^ _ 0.96/Z). 



The last line is a numerical fit to the previous one. 
{Y^pZpUp) « 1.522^1-1/12(1 _ 0.48/Z), so 0.445 = 
0.192 X (1.522)2). In the exchange term: Z^^ is bet- 
ter than 0.96 is close to c in ZS and MS CDE fits. 
The apparently awkward 1-1/12 exponent comes out of 
the fit It will be seen to be natural once we extract b 
from hw-^ calculated as in [Q, Eq. (2-157)], but treating 
separately neutrons and protons: 



35.59(2Z) 



1/3 



■ 0.445 = 



0.860^1/*^ 
''o P 



(rl) — 3rlp'^/5 to bring Eq. 



(7) 



6^ to 



where we have used 

the form of Eq. @). The factor in Z^/^ = (Z2)V 
(Z(Z — l))i/i2 in (;itj)i/2 conveniently corrects the "awk- 
ward exponent" and both equations become identical to 
within a 0.5% discrepancy in the d coefficient: 0.860/ro 
in Eq. (|) vs. 0.864/ro in Eq. Therefore, Ecm is 
essentially exact to within the averaging of L^. A full 
treatment of this term is of obvious interest. 
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Let us turn to the CED in the pf shell. They are 
differences in expectation values of Vc for excitation en- 
ergies; CEDj = E'^{Z,N,J) - E^{N,Z,J) {Z > N). 
The wavefunctions are obtained through standard shell 
model calculations Jl5|,p|-[7| (they depend very little on 
the interaction: KB3, KB3G or FPD6). First we sepa- 
rate monopole and multipole pieces: 



CEDj = ^{Vc,n)j + ^{Vcm)j 



(8) 



(ycm)j is proportional to the difference of (inverse) radii 
between a J yrast and the ground state: — R^^ w 
{Rj — Rq)/Rq). Since Rj is very nearly the same for both 
members of the mirror pair (remember: the neutron skin 
is small), it will be proportional to the average neutron 
plus proton occupancies for the individual orbits, which 
we denote by {mk)j/'2, with nik = + Uk {uk is the 
number of neutrons). Now: to good approximation, the 
Rq denominator is a constant over the region of interest 
{A = 47-51) , furthermore, it is reasonable to assume that 
radii of the pf orbits depend only on I, and, finally, the 
Pi/2 occupancy is always negligible. Therefore, {Vcm)j 
can be taken to depend only on lyvn^^i^ j 12. 

The multipole contribution A(Vc'm),/ is given by the 
expectation values of the effective Coulomb potential in 
the p] shell. Then 



CED,7 = a™(mp3,,)j + A(F^^;)j 



(9) 



The value of be estimated by noting that the sin- 

gle particle P3/2 state in "'^Sc is 200 keV below its ana- 
logue in ^^Ca. This number comes from two effects: a 
larger radius that depresses the P3/2 orbit, and the single 
particle I? term in Fig. |l|that depresses the /7/2 ground 
state orbit. The latter effect is readily found to be ~ 150 
keV by expanding Up + .OIL^)^ around the A ^ AO 
closed shell and using the numbers in, and after Eq. (|^). 
Then, a™ « (.200 + .150)/2 = 0.175 MeV. Note that the 
single particle contribution in is proportional to the 
difference of proton and neutron occupancies. It is im- 
portant in A = 41, but typically ten times smaller than 
the radial effect in A =47-51, so we have neglected it. 

The available information on A(Vppj)j involves only 
the f'j/2 matrix elements extracted from the "^^Ti-^^Ca 



pair, which yields (7 = /y/a) V2:Af7777 = ^ch, 
116.9, 10.9, -59.1) keV, for J= 0, 2, 4, 6 respectively, 
to be compared with the ho values yc%/-,~ (81.6, 24.6, 
-6.4, -11.4) keV. Since we are dealing with Vcm, the cen- 
troids V77 (Eq. (^)) have been removed. They are very 
close (304 keV for A=42, 308 keV for ho) because Ecm 
depends on conserved quantities that cannot be renormal- 
ized. The multipole matrix elements are very different, 
and the data unequivocally prefer the A = 42 set [^1). 
The strategy adopted in these references was to use a 
^Cpf ■with ^cfy2 ' l^^^Pii^S ^Cpf the other matrix el- 
ements, which turned out to be almost irrelevant: the 
Vni^ set by itself accounts for the full Vf^'^l chosen in 
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this way. The results alternated between agreement and 
distorsion of the observed patterns. 

As there is no justification in accepting an enormous 
renormalization for the /y/a orbits and leave the rest of 
the interaction unchanged, we shall attempt a more gen- 
eral treatment, by exploring the possibility of writing an 
effective interaction solely in terms of V^'f^^ — properly 

renormalized to give a plausible account of the full V'^p^ ■ 
First we check that the program can be enacted for the 
ho set. We try 



A{vS;f)j = bA{VSl^Jj + a< 



(10) 
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FIG. 2. Calculated CED normalizations 

Eq. ( |lO| ) is only a numerical recipe, but it works well, 
as seen in Fig. |[ representative of the quality of the ad- 
justment in the four cases (parameters in Tab |). Next, 
assume that the prescription applies to the renormalized 
case and try to invent an effective interaction consistent 
with the data and respecting the condition that the a„i 
parameter in Eq. be a constant that should be deter- 
mined carefully: There is no room for invention here. 

TABLE I. Adopted parameters in Eq. (Q and Eq. (^ 
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by a factor (41/50)2/3 « 0.88 to account for the de- 
nominator (after Eq. (H)). Therefore, we set Um = 0.15. 
Now choose 6m = b/2, aM = 1.8a. Eq. (|n) with these 
parameters (Tab ||) yields CED that in Fig^ are seen too 
agree well, even very well, with experiment. The mild ex- 
ception is A = 50, discussed in detail in Ref. |7|, which 
contains a heuristic introduction to our CED results . 

The monopole Vcm and niultipole Vcm contributions, 
shown separately in Fig.^, indicate that the latter re- 
produces only roughly the experimental patterns: The 
addition of Vcm is indispensable to bring quantitative 
agreement. It is especially worth noting that the strong 
signature effect in the A — 49 band is erased in the CED 
by the out-of-phase Vcm and Vcm- Conversely, the sig- 
nature staggering is enhanced in ^ = 51. 

The monopole contribution provides valuable informa- 
tion about the evolution of yrast radii. As a consequence, 
the use — and even the validity — of the schematic multi- 
pole term (the "invention" ) must be assessed by the focus 
it brings to the monopole one, which must be present in 
a form very close to that in Eq. (^). 

To conclude: once the NS anomaly is resolved, the 
Coulomb field fulfills the — long held — hope of providing 
information about radii not directly measured. The 
terms offers intriguing prospects. A complete analysis of 
Vcm, including non-diagonal contributions, is in order to 
estimate isospin impurities. The renormalization of Vcm 
remains an open problem. 

This work owes much to a stay of AZ at the UAM, 
made possible by a scholarship of the BBVA foundation. 
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FIG. 3. Experimental fl-H and calculated CED for the 
pairs 4^Cr-4^V, ^«Mn-^«Cr, ^^Fe-S^Cr, and ^^Fe-^^Mn. 



Let us write 
CEDj = Vctti + Vcm 



(11) 



For A « 50, 



i("^p3/2)J + bM A{V^y^^Jj + aMimp^^Jj 
•,m — estimated at A = 41 — must be reduced 
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